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Abstract 

Minimum phase functions are fundamental in a range of applica- 
tions, including control theory, communication theory and signal pro- 
cessing. A basic mathematical challenge that arises in the context of 
geophysical imaging is to understand the structure of linear operators 
preserving the class of minimum phase functions. The heart of the 
matter is an inverse problem: to reconstruct an unknown minimum 
phase preserving operator from its value on a limited set of test func- 
tions. This entails, as a preliminary step, ascertaining sets of test 
functions that determine the operator, as well as the derivation of a 
corresponding reconstruction scheme. In the present paper we exploit 
a recent breakthrough in the theory of stable polynomials to solve the 
stated inverse problem completely. We prove that a minimum phase 
preserving operator on the half line can be reconstructed from data 
consisting of its value on precisely two test functions. And we derive 
an explicit integral representation of the unknown operator in terms 
of this data. A remarkable corollary of the solution is that if a linear 
minimum phase preserving operator has rank at least two, then it is 
necessarily injective. 
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1 Introduction 

In the space L^(M+), which is the basic model for causal signals of a con- 
tinuous variable, the term minimum phase refers to functions / which are 
determined by their power spectrum, i.e. the modulus of the Fourier 
transform, and which among all functions having a given power spectrum 
maximize the partial energies Jq for all T > 0. Minimum phase func- 
tions are fundamental in myriad contexts, including control theory, commu- 
nication theory and signal processing. This paper presents a full solution to 
the following inverse problem, which is motivated by geophysical consider- 
ations. Let A be an unknown linear operator on L^(M+) with the property 
that the function Af is the translate of a minimum phase function when- 
ever / is the translate of a minimum phase function. Which sets of functions 
S C L^(M4.) determine Al And how can A be reconstructed from its values 
AS on such sets? The main result of the paper. Theorem [6] in Section HI 
gives an explicit formulation for such an operator A in terms of its value at 
the two functions, e~*(l — t) and e~*t. The proof uses a new result in the 
theory of stable polynomials that characterizes all linear operators mapping 
polynomials having zeros outside the open unit disk to zero-free analytic 
functions on the disk. Hardy spaces on the half plane and disk serve as the 
basic framework for our arguments. 

In geophysics there is a long standing idea — see [7], for example — that 
in a horizontally stratified absorptive earth with vertically traveling plane 
compressional waves (represented as functions of travel time), transmitted 
waves are translates of minimum phase functions. And there have been field 
experiments, [8], supporting the assertion that, for instance, the source sig- 
nature of a dynamite blast is minimum phase. This provides the conceptual 
basis for the inverse problem considered in the present article, supposing 
that the aforementioned minimum phase waveforms are the output of a lin- 
ear operator that appropriately encodes the material properties of the earth. 
In the present article no attempt is made to analyze the physical model of a 
layered earth, nor to explicitly represent solution operators for the attendant 
governing equations. Rather, the minimum phase hypothesis is taken as a 
starting point, and the aim is to investigate its consequences. However, the 
precise results obtained in the present paper offer the prospect of practical 
significance for geophysical imaging, as well as implications for the wider 
arena in which minimum phase signals play an important role. 

The paper is organized as follows. Section [21 below, assembles the nec- 
essary notation and background material required to solve the main inverse 
problem. The background material is divided into four subsections: Sec- 
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tion 12.11 covers the notation and standard factorization theorems for the 
Hardy-Hilbert spaces on the half plane and disk; Section 12.21 establishes an 
isometric isomorphism between L^(R_|_) and the Hardy space on the disk, 
which helps to streamline later definitions and proofs; Section 12.31 covers 
the technicalities of minimum phase, as it relates to both the discrete and 
continuous settings; and Section 12.41 quotes two results from the theory of 
stable polynomials that are needed to prove the main results. 

There are two main steps toward solving the main inverse problem. The 
first step is to derive a structural characterization of linear operators on 
L2(M+) that preserve translated minimum phase operators. This is carried 
out in Section [3l where Theorem [3] describes minimum phase preserving 
operators in terms of product-composition operators on Hardy space on 
the disk, and Theorem H] describes them in terms of product-composition 
operators on Hardy space on the half plane. 

The final step toward solving the main inverse problem is to reconstruct 
an unknown minimum phase operator in terms of test data. This is done 
in Section HI where two explicit ways of identifying an unknown minimum 
phase preserving operator are given, in Theorems [5] and [6] respectively. The 
characterization from the previous section is incorporated into the statement 
of Theorem [6l making the latter the main result of the paper. 

In Section m the case of untranslated minimum phase signals is handled 
briefly in Section 15.11 Then Section 15.21 points out a remarkable corollary 
of the main result, that a linear operator that preserves minimum phase is 
necessarily injective, with the exception of the rank 1 case. 

Lastly, Section [6] briefly concludes the paper. 

2 Notation and background results 

2.1 The Hardy-Hilbert spaces on the half plane and disk 

Let H^_^ denote the Hardy space consisting of all analytic functions on 
the open right half plane C+ whose square integrals along lines of constant 
positive real part are uniformly bounded, and let denote the Hardy space 
of all analytic functions on the open unit disk B whose Taylor coefficients are 
square summable. Both and are separable Hilbert spaces. Indeed 
the operator 



defined by 
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for every F £ ^^'^ z G B is a surjective isometry, with inverse 

1 „/l-w' 



(cI>-iG)M = -=— — -G 



^/^T{l + w) \1 + w 

See [3 Chapter 8]. 

A function G E is defined to be outer if G is not the zero function 
and there exists a scalar A of modulus 1 such that for every z G B 

G is defined to be inner if |G(e*^)| = 1 for almost every point e*^ belonging 
to the unit circle S^. It is a fundamental fact that, apart from the zero 
function, every G G H"^ has a factorization of the form 

G — Ginncr Gouter ) 

where Gjnner S H'^ is inner and Gouter £ H"^ is outer; this factorization is 
unique up to multiplication by a scalar in . 

A non-zero function F G -f^c+ defined to be outer if there exists a 
scalar A G 5"*^ such that for every w G C+, 



FH = Aexpfi r y^\og\F{iy)\ 



dy 



1 + r 



Using the fact that 1 + z G -ffj^ is outer, direct calculation shows that a 
function F G H^_^ is outer if and only if its image ^F G is outer. 

The inner-outer factorization in translates via ^ to as follows. 
Given F G let G = GinnerGouter = ^F, and let m : CU{oo} CU{oo} 
denote the Mobius transformation 

m{z) = m~'^{z) = ^—^, (1) 

which maps B onto C+ and vice versa. (In signal processing m{z~^) is 
known as the bilinear transform — see [UJ Chapter 7].) Observe that 

•I" ^(GinnorGoutor) = (Ginncr ° W-) ■ ("I" ^Goutor)- 

This induces the factorization F = Finner-^outer, where the factor 

-^'inner — Ginner ^ 'Ti 



9 September 2011 



PC Gibson, MP Lamoureux 



Operators on L^(M+) 



5 



has modulus 1 on the imaginary line w = iy, and the factor 

-^outer — ^ G'outer 

is outer. Note that the factorization on differs from the one on in 
that the factor Finnerj which is not square integrable on the imaginary line, 
does not belong to H^_^ . 

2.2 The isometry H mapping L^(]R+) onto 
The Fourier-Laplace transform 

C : L\R+) ^ Hl^ 

defined by 

for every / G L^(M+) and w G C+ is a surjective isometry. Its boundary 
function is the Fourier transform of /, denoted Tf; that is, 

{Cf){iy) = {:Ff){y) fora.e. yGM. 

Let T-L : L^(R_|_) — > denote the composition 

H = 

Note that H is a surjective isometry, since each of its factors is. The map 
T-L is given explicitly by the formula 

for every / G L^(M+) and z G B. Its inverse mapping 

H-^ : Hi ^ l2(M+) 

is given by the formula 

for every t G M+, provided the boundary function of is in L^(M); 

otherwise T-l~^G may be computed as an oscillatory integral as per the 
standard theory for Fourier transforms on L^(R). 
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For each nonnegative integer n, let /i„ e denote the nth moment 
function 

Since the functions /in comprise an orthonormal basis for H^, and 'H~^ is 
an isometry, it follows that the functions T-l~^fin form an orthonormal basis 
for L^(M_|_). These may be computed explicitly using the above formula for 
H"^, which, after some work, yields 

(^-Vn)(i) = (-l)"V2e-*L„(2i), 

for each n > 0, where L„ is the nth Laguerre polynomial defined as 

t m 

n!dt"^ ' 

In particular, 

(^-Vo)(i) = V2e-* and (^"VOW = \/2e-*(2i - 1). 

Let 

denote the map defined by the formula 

0/ = {a4~ 0, 

where for each n > 0, 

/"oo poo 

an= f{t){H-^finmdt ={-irV2 f{t)Ln{2t)e-'dt. 
Jo Jo 

Note that, since the sequence of Fourier coefficients {cn}^=Q = FG of a 
function G G is given by the formula 

-I /'27r poo 

cn = (G, Mn) = 7^ / G(e*^)e--^ d0 = / (H-^G)(t)(7{- Vn)(t) dt, 

the map d simply associates the sequence of Fourier coefficients {onj^o ~ 
FT-Lf to a given continuous signal /. Since the transform 

is an isometry, Z) is therefore an isometry too. 
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2.3 Minimum phase functions 

The hiverse of the transform F : — )• is the z-transform 



defined by the formula 



(Za) (z) = ^a„z" 



n=0 



for every sequence a = {a„}^Q G ^^(^+) and z G ID). A causal discrete 
signal in is defined to be minimum phase if its z-transform is an outer 

function, or in other words, if the given discrete signal is the sequence of 
Fourier coefficients of an outer function. A minimum phase signal {a„}JJ^Q 
has maximally front loaded energy, in the sense that if {6„}5JLo ^ ^^(^+) 
satisfies 



E 



for every z G 



n=0 

D, then for every > 0, 
N 



oo 

E 

n=0 



N 



CLnZ 



n=0 



n=0 



See [6]. 

A function / G L'^(M_|_) is defined to be minimum phase if its image 
Hf G is outer. As in the discrete case, if / is minimum phase then its 
energy is maximally front loaded as follows. For every g G L^(M_|_) such that 
\Tg\ = l-^/l, and for every T > 0, 

\f{t)\^dt> r \g{t)\' dt. 
Jo 

In the continuous setting the connection between minimum phase and the 
energy being maximally front loaded is much more delicate than in the 
discrete setting, where the two properties are well known to be equivalent. 
Non-uniqueness of solutions to the Pauli problem — see[l] — bedevils the con- 
tinuous case; this will be treated in a separate paper. 

The symbol Ai denotes the set of all minimum phase functions in L^(M-i-). 
Given / G and r > 0, the translate 



{Trfm 



f{t-r) 




if t > T 
if t < T 



9 September 2011 



PC Gibson, MP Lamoureux 



Operators on -L^(M+) 



8 



is minimum phase only if r = 0. Nevertheless, from the signal processing 
point of view, translates of minimum phase functions are of similar interest 
to minimum phase signals; they are also front loaded, but with respect to 
a given delay r > 0. The set of all translates of minimum phase signals is 
denoted Ai. That is, 

M = {T^f I / € X and r > 0} . 
The foregoing definitions of minimum phase imply the following. 
Theorem 1 As in Section [2^ let 



be defined by the formula 
where for each n> 0, 



5/ = {anr 



n=0' 



an = / f{t){n-^tinm dt = i-irV2 / /(t)L„(2f)e-* dt. 
Jo Jo 

Then d is a surjective isometry with the property that t)/ is minimum phase 
if and only if f is minimum phase. 

Proof. As noted at the end of Section [221 f = FT-L is an isometry since each 
of 7i and F is. By definition, a function / G L^(M-|_) is minimum phase if 
and only if Tif is an outer function, which is true if and only if its sequence 
of Fourier coefficients 

f / = FTif 

is minimum phase. I 



2.4 Stable polynomials and stable preserving operators 

The arguments given in the present paper use results from [2j concerning 
stable polynomials, necessitating some notation as follows. Given a set 
0, C C, the algebra of all analytic functions / : J7 — )• C is denoted by 
A{0,). Given sets O C fi' C C, a function / : fi' — )• C is said to be 
stable if f{z) ^ for every z £ 0,. The set of all 0-stable functions in 
A{0,') is denoted by S{^1). The symbol V{il.) denotes the set of all il-stable 
polynomials. Given tp G A{0,), 

: A{n) A{n) 
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denotes the operation of multiplication by -0, 

iM^f)iz)=^Piz)f{z) 
for every / S A{i^) and z £ Q. Given an analytic function (p : Q ^ Q, 

: Ain) Ain) 

denotes the operation of (right) composition with ip, 

{C^f){z) = fo^{z) = f{ip{z)) 

for every / G A{^) and z £ VL. The symbol C[z\ denotes the set of all 
univariate polynomials over C. 

Theorem 2 (from [2^, Theorem 2]) LetVL he a non-empty connected 
open set. If a linear map A : Clz] — )• A{p) has the property that 

A{V{I})) C S{a) U {0}, 

then either: 

1. there exist a function ip € S{Q) and a linear functional : C[z] ^ C 
such that A{f) = iy{f)'ilj, for all f € C[z]; or 

2. there exist a function ■0 G S{Q) and a non- constant function ip G ^(D), 
where (p{^) C B, such that A = M^C^. 

Note that there is no topology attributed to the vector space C[z], and 
hence no hypothesis of continuity in the above theorem; the maps in question 
are simply linear, with no additional restrictions. By contrast, the next 
result requires boundedness. Let O denote the set of all outer functions in 

Lemma 1 (from [2l Lemma 4]) ^4 bounded linear functional u : — )• C 
satisfies 

v{0) C C\{0} 

if and only if there exist a point zq GO and a scalar a G C \ {0} such that 
for all f £ H^, 

v{f) = af{zo). 
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3 The characterization problem 

Theorem [2] makes it possible to give describe explicitly the structure of a 
linear operator A : L^(M+) -^^^(1^+) for which A{M) C M. 

Theorem 3 Let A : L^(M-(_) — t- L^(IR+) be a bounded linear operator such 
that A{M^) C M.. Then there exist a function ip G 5(D)ni7^ and an analytic 
function 93 : D — )• ID such that 

A = n-^M^c^n. 

Proof. Consider the associated operator 

B = UAW^ : Hi Hi. 

Noting that every stable polynomial p G V(Ji) is an outer function, the 
hypothesis that A preserves translates of minimum phase functions implies 
that 

AU-^p G M. 

Letting g £ M have the form T^f, where r > and f £ A4, the definition 
of Ti yields that 

{ng)iz)=exp(T^) {nf){z). 



Since / is minimum phase, the function Hf is outer. But if r > 0, then Hg 
is not outer, since the function e ^+1 is singular inner; see [5, Section 2.6]. 

^ z-l „ 

Nevertheless, e ^+1 is B-stable, as is every outer function in Hf^, and there- 
fore Hg £ 5(B). This shows that 

B{V{B)) C cS(B), 

so that B satisfies the hypothesis of Theorem [2] with i7 = D. Thus B is 
either a product-composition operator or a rank one operator. 

By hypothesis B does not map any outer function to 0. This implies that 
if i? is a rank one operator of the form given in the first part of Theorem [2l 
then the associated linear functional v satisfies the hypothesis of Lemma [TJ 
The conclusion of Lemma [1] then implies that B has the form 

B = M^Cip, 

with ip{z) = zq a constant function. This shows that B may be represented 
as a product-composition operator whatever its rank, establishing the de- 
sired characterization of the original operator A. I 
At this point the issue of sufficiency naturally arises, and there are two 
basic questions. What conditions on ip and (p guarantee that: 
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1. M^C^ : — )• is a bounded operator; 

2. the associated operator A = T-L^^ M^C^V. maps A4 into itself? 

Regarding question 1, for any analytic function 93 : D — )■ D and any ip S H^, 
the product-composition operator M^C^p is a bounded linear operator. See 
[5l Theorem 5.1.5]. But these conditions are not necessary in that ip need not 
be bounded. See jH Section 4] for a more detailed discussion. Concerning 
question 2, the operator A = H^^M^C^pH preserves M. if the product- 
composition operator M^dp preserves T-L{A4.), which consists of functions of 
the form 

where p G is outer and r > 0. If ^/^ G 'H{M) and either Lp{z) = z 01 the 
function 

1 + 

has an analytic extension to the closed disk D, then M^C^ is guaranteed 
to preserve T-L{A4). These conditions are not particularly stringent, and the 
corresponding class of operators is large. 

Returning to the earlier thread, A = H^^M^C^pT-L has the integral rep- 
resentation 

{Af){t) = - / / r,{y)e'^y^-''') fir) drdy, 



where 



and 




The above formulation is somewhat awkward on account of its representation 
in terms of functions on the unit disk. An alternative is to express A = 
T-L~^ M^C^T-L in terms of functions on the right half plane, as follows. 
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Theorem 4 Let A : — t- L^(]R_^) 6e a continuous linear operator 
such that that A{Ai) C A4. Then there exist a function a G 5(C4.) n H^_^ 
and an analytic function ^ : C+ — t- C+ such that 

A = C-^M^C^C, 

with K = \/27r(l + ^)a. 

Proof. By Theorem [3] the operator A has the form 

where ^ G 5(ID)) n i^^ and : D — 7- P is analytic. The operator $ ^MipCip^^ 
is itself a product-composition operator, as follows. For any F G H^_^ , 

2^^ ^/1-V9(z) 



l + (^(z), 
i*" o m o 

-v/7r'0(^)(l + m o (p)(F o m o (p){z). 



2 

1 + 



Therefore 

+ m o o m)F o m o (f o m 

= M^C^F, 

where 

Since (/? : D — t- D is analytic, it follows that ^ : C-(- — t- C+ is analytic too. 
Note that carries B-stable functions in to C+-stable functions in 
H^^, and set 

so that a G -fT^.^ is C+-stable, just as tp is B-stable, the factor of l/\/2 being 
for later convenience. In terms of this notation, 

A = C'^M^C^C where k = V2^{1 + 0^^, 
completing the proof. I 
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The formulation of Theorem U] is preferable to that of Theorem El since 
an operator A = C~^Mf^C^C with k = -v/27r(l + has an integral repre- 
sentation of the comparatively simple form 

1 i^oo poo 

{Af){t) = ^ / / a{iy){l + C(iy))e--«(*^)+^^V(r) drdy. (2) 

V^Vr Jy=-oo Jt = 

The reason for expressing k = ^/2tt{1 + 0'^ factored form is that, from 
the point of view of test data, a is the primary object. In the next section 
it will be seen that C~^a can be measured directly by evaluation of A at 
the test function e~*, whereas k, must be computed from a combination of 
measurements. 

4 The identification problem 

With Theorem[3]in hand, the problem of identifying a given minimum phase 
preserving operator using test functions becomes straightforward. 

Theorem 5 Let A : L'^{R+) L'^i^+) have the property that A{M) C M. 
Then A can be reconstructed from its values at po{t) = \/2e~* and pi{t) = 
\/2e^*(2t — 1) by means of the formula A = %~^M^Cip%, where 

ilj = nApQ and if = {UApi) / {UApo). 

Proof. Recall that 

po(i) = (-H-Vo)(t) and = (7^-Vi)(t), 

so that 

Apo=^"V and Api=n~^{iljtp). 

I 

Corollary 2 The function tp in the conclusion of Theorem has the form 

friz - 

i^{z) = expl \ i^oiz), 

where r > and ipo € -ffp is outer. 

Proof. The test function po is minimum phase, so by hypothesis ApQ is 
the translate by some r > of a minimum phase function. It follows that 
TiApo has the stated form, just as in the proof of Theorem [3l I 
Note that the functions po and pi in Theorem [5] can be replaced by any 
pair having the same span. The following choice is geared to fit Theorem S] 
and the formulation ([2]), resulting in simple reconstruction formulas. 
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Theorem 6 Every continuous linear operator A : L^(M+) — t- L^(M+) with 
the property that A{Ai) C M is determined by its values Aao and Aai, 
where 

aQ{t) = e~*(l - t) and ai{t) = e~*t 
for all t G M+. Setting a = CAgq + CAai and ^ = LAo^j LAa\, 

A = C'^M^C^C with K = V2^{1 + ^)a; 

equivalently, A has the integral representaion 

iAf){t) = ^ / / a{iy){l + e(iy))e-^«(^^)+*^V(T) drdy. 



~oo J r=0 

Moreover, a necessarily has the form a{w) = e~^'^aQ{w) for some scalar 
e > and some outer function oq G ^c_|_ ' ^''^'^ ^ • ~^ ^+ analytic. 

Proof. Note that 



w 



27r(l + w)2' 
1 



(£cro)(w;) = 

Using these formulas and Theorem [H a straightforward calculation gives the 
desired formulas for a and ^. Since gq+cfi is minimum phase, the hypothesis 
that A preserves translates of minimum phase functions implies that a is 
the Fourier-Laplace transform of a translate by some e > of a minimum 
phase function. It follows that a{w) = e~^"'ao(tt;) for some outer function 

There is no way to identify an unknown operator of the form A = 
M^CipTi. from a single value Af, no matter which test function / G 
L^(M+) is chosen. In this sense Theorems [5] and [6] give results that are best 
possible. 



5 Remarks 

5.1 Operators preserving Ai 

The identification of bounded linear operators A : L^(M+) — )■ L^(M_|_) pre- 
serving minimum phase functions, but not necessarily their translates, can 
be handled using exactly the same methods as above. Indeed without con- 
sidering translates, the arguments become simpler, and the analogue of The- 
orem [6] that results is only slightly different, as follows. 
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Theorem 7 Every continuous linear operator A : L^(M+) — t- L^(M+) with 
the property that A{A4) d M is determined by its values Actq and Aai, 
where 

aQ{t) = e"*(l - t) and ai{t) = e^^t 
for all t G Setting a = CAa^ + llAa\ and = CAa^j LAa\, 

A = C'^M^C^C with K = V2^{1 + ^)a; 

equivalently, A has the integral representaion 



Moreover, a £ H^_^ is necessarily an outer function, and ^ : C+ — t- C+ is 
analytic. 

Thus the only change is that the function a in the integration kernel is 
necessarily outer, whereas before it was allowed to be a modulated outer 
function. Another way to approach the case A{M.) C is to use the 
minimum phase preserving isometry 



defined in Section 12.21 and appeal to the corresponding theorem for the 
discrete case, |2i Theorem 4]. The mapping D does not carry translates of 
continuous signals to translates of discrete signals, so by contrast Theorem [6] 
cannot be proved using its discrete analogue. Indeed, whereas Theorem [6] 
needs the full force of Theorem [2 its discrete analogue was first proved in 
[H Theorem 4] using less powerful methods. 

5.2 The hypothesis of boundedness 

It is a remarkable fact that if the operator A : L^(M+) — )• L'^(M.^) is assumed 
to have rank at least two, then the hypothesis of boundedness can be dropped 
from Theorems O HI O [6] and [71 Boundedness was only used when invoking 
Lemma [11 which applies to the rank 1 case. Moreover, in the case of rank 
at least two, the property that A{M) C M or A{M) C M forces A to 
be injective, a property implied by the product-composition structure of 




UAU-^. See Corollary 2]. 
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6 Conclusion 

This paper's main results, Theorems [6] and [7] in Sections H] and 15.11 solve 
completely the natural inverse problem concerning linear operators on the 
half line that preserve the class of minimum phase — or translated mini- 
mum phase — signals. They show that every such operator is conjugate via 
the Fourier-Laplace transform to a product-composition operator on Hardy 
space, that precisely two test functions are required to reconstruct such an 
operator, and that there is a simple reconstruction formula in terms of the 
exponentially damped linear test functions e~*(l — t) and e~*t. Recent de- 
velopments in the theory of stable polynomials played a crucial role in the 
derivation of these results. 

The inverse problem of identifying a minimum phase preserving oper- 
ator from its output was motivated by considerations in seismic imaging, 
where seismic traces are widely assumed to be minimum phase functions, or 
translates of minimum phase functions. Results in the present paper show 
that if the minimum phase hypothesis is valid for seismic data, then cur- 
rent techniques such as Wiener deconvolution are inadequate in the most 
general setting. A non-stationary operator representing the process of an 
input wavelet travelling through an interval of rock, if it preserves minimum 
phase, cannot be determined from a single test function. And therefore 
an input wavelet cannot be computed from a single seismic trace — even 
with non-stationary methods; at least two traces, arising from two different 
source functions are required to determine the operator, and thus to invert 
it. This offers the exciting possibility that a re-evaluation of the seismic 
deconvolution techniques in light of the Theorem [6] could lead to substantial 
improvement. In any case the very precise results in the present paper shine 
a spotlight on the minimum phase hypothesis itself, and call for a better 
understanding of its origins in models of layered media. 
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